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Abstract 

The astrometric signature imposed by a planet on its primary increases substantially 
towards longer periods (oc P 2 / 3 ), so that long-period planets can be more easily detected, 
in principle. For example, a one Solar-mass (M ) star would be pulled by roughly 1 milli- 
arcsec (mas) by a one Jupiter-mass (Mj) planet with a period of one-hundred years at a 
distance of 20 pc [cf. eqn. ([3]) below]. Such position accuracies can now be obtained with 
both ground-based and space-based telescopes. The difficulty was that it often takes many 
decades before a detectable position shift will occ ur. However, by the time the next gener- 
ation of astrometric missions such as SIA^ [e.g., ( Edberg et al. 2005 )] will be taking data , 
several decades will have past since the first astrometric mission, HIPPARCOS ( 1ESA 1119971 ). 

Here we propose to use a new astrometric method that employs a future, highly accurate 
SIM Quick-Look (SQL) survey and HIPPARCOS data taken twenty years prior. Using a 
conservative position error for SIM of 10 /ias, this method enables the detection and char- 
acterization of "Solar-system analogs" (SOSAs) with periods up to 80 (165) years for 1 (10) 
M j companions. Employing the standard SIM error of 4 //as, this period range is extended 
by a factor of two to four. We might expect the PDF to turn over in this period regime. 

Because many tens of thousands nearby stars can be surveyed this way for a modest ex- 
penditure of SIM time and SOSAs may be quite abundant, we expect to find many hundreds 
of extra-solar planets with long-period orbits. Such a data set would nicely complement 
the short-period systems found by the radial-velocity method. Brown dwarfs and low-mass 
stellar companions can be found and characterized if their periods are shorter than about 
500 years. This data set will provide invaluable constraints on models of planet formation, as 
well as a database for systems where the location of the giant planets allow for the formation 
of low-mass planets in the habitable zone. 

1. Introduction 

Thi s white paper is ba sed on very recent work summarized in a r eview paper on SIM- 
science ( jUnwin et al. 1 120071 ). while we present more details elsewhere tolling 1 120071 ). 

Our current knowledge of the demography of extra-solar planetary systems is mostly a 
result of long-term radial velocity (RV) monitoring of nearby (mostly) FGK main-sequence 
(MS) stars. The results are spectacular with over 200 suspected planets in >171 systems 
( {Schneider 1120061 ). with most of the planets in short-period orbits. Only 10% of the observed 
planets have periods exceeding 5 years, while just one (0.5%) has a period slightly longer 
than the period of Jupiter (11.9 yr). We use the P DF of extra-solar giant planets (ESGPs) of 



Tabachnik fc Tremaine I (12002. hereafter TT2002I ) but scaled-u p by a factor 1.6 to account 



for current understanding the ESGP frequency [e.g., ISozzetti I (120051 )]. The updated PDF 
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indicates that the period range between 5 years and the maximum currently known period 
should account for 25% of the total number of ESGPs rather than the observed 10%. In 
fact, the PDF of ESGPs increases towards longer periods (P) so that systems dominated 
by long-period planets such as the Solar system may be quite common, and we will use it 
to estimate the frequency of solar system analogs (SOSAs). However, the PDF^sgp has to 
turn over at some period to yield a finite integrated probability If we had to guess where 
the planetary PDF might turn over, we might pick th e period where the stellar PDF turns 
over, or about 170 years (jDuquennoy k. Mayor Ill99ll ). 

We loosely define a SOSA as a system with a (single) planet in the mass range between 
Jupiter and Uranus/Neptune (~0.05 Mj) and with periods between 11.9 years (Pj up iter) 
and 165 years (PNeptune)- Integrating the PDF ESGP over these ranges, we find that 12.6% 
of systems would be solar-system analogs. If we consider the group of long-period planets 
that can be detected astrometrically, we need to consider more massive systems with masses 
(M) between 1 and 13 Mj. We call such systems heavy SOSAs, or HOSAs. The PDF esgp 
predicts that such systems make up about 20% of the total number of planetary systems 
with periods up to 165 years, and occur around 7.9% of apparently single stars^l 

NASA's SIM Planet Quest can det e ct ext ra-solar planets weighing several times the mass 



of the Earth [e.g., Catanzarite et al."l ( 20061)]. If m u ltiple planets exist, their properties can 



also be determined with SIM [e.g., ISozzetti et al.l (120031 ); iFord I (120061 )]. The mission-end 



astrometric accuracy of ESA's GAIA astrometric mission is about twenty times worse than 
SIMs, rendering it not very useful for the project described here. 

2. Finding Solar-System Analogs 

There are several astrometric methods that can be used to identify "long-period" com- 
panions of stars. These methods are based on the fact that the actual proper motion (/x) is 
non-linear if it contains a contribution from the reflex motion of the primary being orbited by 
a companion. In the "A/z method" a substantial difference between /i values from a "short- 
term" catalog such as HIPPARCOS and t hose from a "long-term" p r oper motion catalog 



Makarov fc Kaplan 



such as TYCHO-2 is indicative of binar ity dWielen et al. II1999I . l2000h . 
( 12005; hereafter referred to as MK2005I ) find that both the period and mass can be esti- 
mated fro m the acceleration (jx) and t he jerk (/i), but only z/the "long-term" proper motion 
is known. Kaplan fc Makarov I (120031 ) developed a method that is appropriate for objects 
with periods up to twice the mission duration (i.e., up to 10 to 20 years for SIM). 



2 Recent work [Cumming et al. (in preparation) as previewed by Butler et al. ( 2006h ] indicates that the 
mass function declines more rapidly towards higher masses: dN/dM oc M -19 rather than dN/dM oc M~ 1A 
as derived by TT2002. As a result, the number and relative frequency of HOSAS would decrease, while their 
detectability is unchanged. Because we deal with detectability, we will use the re-scaled TT2002 PDF. 
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2.1. Past and Future Astrometry 

Here we propose a new method to identify and quantify long-period systems compris- 
ing planetary, brown-dwarf (BD) and main-sequence (MS) companions. Our method uses 
astrometry of earlier epochs but uses the positions rather than the proper motions. The 
other component is a future, highly accurate astrometric mission such as SIM. The idea is 
to fit the SIM data for a given star with a simple astrometric models [e.g., linear, quadratic, 
etc.], and use this model to predict the position of the star at the HIPPARCOS epoch (t#). 
We assume that SIM data will be from 2013.5, leading to and epoch difference of 22 years. 
We assume that the HIPPARCOS data are accurate to ±1 mas. However, this accuracy 
can be improved upon in the future by more careful modeling of the systematic effects, or 
by using the much improved GAIA reference frame to defin e the frame at the HIPPARCOS 
epoch. The former method has been applied recently by Ivan Leeuwen fc Fantino I (120051 ) 
who reduced the HIPPARCOS errors by almost a f actor of three, w hile the latter has been 
used in the construction of the TYCHO-2 catalog (IHog et al.l 120001 ) . 

3. Model Details 

We assume circular, face-on orbits and neglect all the details associated with orbit 
fitting. The work of MK2005 indicates that the results will not be very sensitive to these 
simplification. We also assume that the secondary is "dark," so that the photocenter tracks 
the motion of the primary. The position (z) of the photocenter of the primary is thus a 
function of: 1) the position (zq) at time t — 0, 2) the proper motion of the barycenter ([i z ,b), 
3) the semi-major axis of the orbit of the primary (a D ), 4) the orbital period and phase (f>, 
and 5) the distance (d pc ) in pc. (We use z as shorthand for either x or y). With the period 
in years, the total mass (M to t) in M and the mass of the companion (Mc\j) in Mj, we find: 
x(t) = x t=0 + fi x ,Bt + X (t) y(t) = y t=0 + + Y (t) (1) 
X Q (t) = a Q cos(27rf/P + 0) Y Q (t) = a Q sin (2vrf/P + 0) (2) 
0.9547 / P \ 2/3 „ r 

~ . * /2tt\ 1 /2tt\ 2 2 1 /2tt\ 3 , 

X {t) = a 0tC(P -\ — \ a oM t--\ — 1 a 0tC(P t +- 1 — 1 a D ^t (4) 

Y (t) = a OiS0 + I — I a ,c<f>t- - I -p I a 0)S<j) t - - I — J a Q ^t (5) 

with a 0>C( ^ = a D cos ((f)) and a 0iS( ^ = a D sin ((f)), and where we expand the position change Z a (t) 
due to the orbit to third order to arrive at eqns. (T4]) and ([5]). We identify the orbit-induced 
position (z Q ), the proper motion, acceleration and jerk as the coefficients of the t-terms with 
powers 0,1,2 and 3, respectively. X and Y are a 3 rd -order, orbit-based astrometric model. 
On the other hand, the observed trajectory can be fit by a polynomial up to n th order: 

ZF,SIAi(t) ~ Zo,F,SIM + Z\,F,SIM t + Z 2 ^ S IM r + Z3,F,SIM t S + C(t 4 ) , (6) 
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where the subscript "F, SIM" indicates that the fit is performed to the SIMdaXa only. Zf,sim 
can be evaluated at any previous epoch and compared with the the observed position at that 
epoch. The position error [S z (t)] on zp,siM depends on the accuracy of the fit and strongly 
on the epoch difference. The difference between the true position at the HIPPARCOS epoch 
and the SIM prediction is given by A z {th) = zh — zf,sim( t h), while the significance of 
A z (th) is readily computed. In order to make a significant detection, A z {r H ) has to be 
smaller than the errors on both the SIM prediction and the HIPPARCOS position. 

While zp^siM can fit the space motion during the SIM observing span extremely well, 
the extrapolation of the model to the HIPPARCOS epoch can lead to large A z values when 
the primary has a companion. In general, large A z (th) values indicate heavy companions, 
while small values indicate either no companions at all or a low-mass companion. 

3.1. Period and Mass Estimates from A z [t h ) 
Ideally, one would like to know the motion of the barycenter so that it could be sub- 
tracted from ZpsiM(t) to yield the orbital contribution. In that case, the period follows from 
the ratio of the coefficients of eqns. (jlj) and (jSJ), and would be independent of orbital phase 
and inclination. Unfortunately, because we do not know /j,b, this method can not be used. 

Alternatively, it is possible to eliminate the phase effects is by combining the x and y 
positions differences, at least for face-on circular orbits. Initial investigations indicate that 
the effects of inclination (i) are not all that large, as long as % < 45°. Keeping in mind that 
the results are indicative rather than definite, we proceed w ith circular, face-on orbits. The 
position differences can be found analytically tolling 1120071 ). and read: 



A x = x(t h ) - X(t h ) A y = y{r H ) - Y(r H ) A xy = J A 2 X + A 2 (7) 



2 2 

= J2 I 1 - 2s «P + 2(1 - c«)p 2 } & A%^ = ^ 



4 + c a p 2 - 2s a p 3 + 2(1 - c a )p 4 



with a = 2ttth, s a = sin (a/P), c a = cos(a/P), and p = P/(2ttth), and where the "// + .. ." 
subscripts indicate that an expansion of the orbital motion is used that includes all listed 
components. The position differences can be ratioed to yield a period estimator: 



= kt h ~ 3n „ . (9) 



A xy ^ r P P < T H 

which is accurate for either short or long periods. In the intermediate regime, P oscillates due 
to the trigonometric terms in eqns. (181) . Once the period has been estimated, the companion 
mass follows from solving either of the A xy relations for a Q [and hence mass via eqn. (J3])]. 

The proper motion of the barycenter is not important for this method because it does not 
matter how the observed proper motion is divided between the center-of-mass- and orbital 
components. The SIM model is good because it predicts the observed positions at the SIM 
epoch, while any position difference at the HIPPARCOS epoch depends only on the orbital 
parameters and th, so that A xy can be calculated employing the orbital parameters only. 
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We have performed extensive numerical sim ulations to test analytical relations for the 
position differences derived above (jOUing 1 120071 ). Our modeling of the system comprises an 
implementation of equations with an arbitrary barycentric motion and a periodic signal 
in both coordinates (with random phases). We use this model to predict the position at the 
HIPPARCOS epoch. We then generate, in Monte-Carlo fashion, many random realizations 
of the model which are fitted by a polynomial to the SIM positions only. The so-determined 
SIM astro metric model is extrapolated to the HIPPARCOS epoch to yield A xy . We perform: 
1) a first order fit to compute A xyttl , and 2) a second-order fit for A xy ^ + ^. The numerical 
results are virtually identical to our analytical predictions [eqn. (jHJ)]. 

4. Results 

We ran simulations that might be relevant for the SIM quick-look survey (^5] below): 
5 observations per coordinate per star during a period of 18 months. Conservatively, we 
assume SIM position errors of 10 /xas per observation. The results are presented in Figured] 
where we plot two A xy metrics as determined from the first-order fit (abscissa) and the 
second-order fit (ordinate). In this figure, systems with a given period but with varying mass 
fall along diagonal lines from the lower-left to the upper-right (drawn linesjfl In Figure [IJ 
constant-mass systems are indicated by the dashed lines. The "bunching up" of the lines 
around A xy ^ + ^ = 5 mas is due to the limited accuracy of the SIM-based determination of 
the acceleration, and this limit is used to generate the thick horizontal line (at 3 times this 
level) . Thus the SIM measurements are the limiting factor for the determination of A xy ^ +/:i 
rather than the HIPPARCOS error. This is the reason why the much less accurate GAIA 
data would not be very useful for this application. Figure [1] indicates that the companion's 
mass and period can be determined in a large region of parameter space. The residuals 
suggest that the orbits of a 1 (10) Mj ESGPs can be characterized up to periods of 10 (80) 
years, while this is possible for stellar companions up to 500 yearfl If we use the expected 
SIM accuracy of 4 //as, this period range is extended by a factor 2-4. 

5. A SIM Quick-Look Survey for HOSAs 

Figure [T] also indicates that a few highly accurate observations suffice to identify plan- 
etary, MS and BD companions. Such data could be generated by a SIM quick-look (SQL) 
survey. The targets are bright HIPPARCOS stars, so we assume that an SQL observation 
can be achieved in one minute per position per baseline. Thus, one-thousand stars can be 
done in 10,000 minutes (6.9 days), so that several thousand stars can be included in an SQL 
survey without impacting the overall SIM mission significantly. Given that the PDF for 



3 



The well-known period-mass degeneracy would result if we were to plot the fitted /i instead of Axy.^+p,. 



4 The l si -order A xy ^ values are significant up to 1,000 years at the Hydrogen burning limit and 4,000 



years for a double star with solar-mass components. 
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Position Difference Predicfions 
HIPPARCCS Position SQL Prediction 
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Fig. 1. — For a series of models with various companion masses (vertically oriented numbers 
in the plot in units of Mj) and orbital periods (in years; horizontal numbers) we plot two 
positions difference computed from our model data. The abscissa is the position difference 
(A X y !fJ ), while ordinate is A xy ^ + ^. The 3-sigma observational limits are indicated by the 
thick vertical and horizontal lines. 



ESGPs predicts heavy solar-system analogs around 7.9% of stars, a survey of ~5,000 stars 
may find 400 HOSAs. Such a sample would firmly establish the PDF in the long-period 
regime and indicate how unique the Solar system really is. 

In order to maximize the yield of HOSAs (and some SOSAs), an SQL program nee ds to 



avoid MS and BD multiples. The subset of 73,000 ARIHIP stars flWielen et al. Il200ll ) that 



show no signs of binarity is a good starting point for the target selection of an SQL survey. 

Those systems that do not show signs of binarity in the SQL+HIPPARCOS survey are 
likely to have either sub-stellar companions or stellar companions with very long-periods. 
Those systems warrant further SIM observations. The SQL follow-up survey of those stars 
with suspected sub-stellar companions would be significantly more sensitive than the SQL 
survey. Figures similar to figured] but with employing the SQL follow-up data indicate (not 
shown) that the ESGPs can be detected with masses as low as 0.1 Mj in 10 year orbits. 
Period estimation for 1 [10] Mj is extended by a factor four [two] (to 40 [160] years). 

6. Conclusions 

A judicial combination of HIPPARCOS data, a SIM quick-look survey and follow-up 
SIM observations at full accuracy can uncover several hundred extra-solar planetary systems 
with periods comparable to the gas giants of the Solar system. Such a program is only 
possible with SIM-like accuracies, and the results would nicely complement radial velocity 
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and imaging surveys. 

Given the importan ce of accurate pre-SIM astrometry, it is sad to realize that the 



canceled FAME mission (jJohnstonll2003l ) would have provided an excellent reference catalog 
for detection and characterization of solar-system analogs. Likewise, it is of pre-eminent 
importance to continue all-sky astrometric programs at intervals of ten to twenty years to 
probe the long-period regime. The required accuracy depends on the desired period- and 
mass ranges, but a survey with an accuracy at the HIP PARC OS level (one-half to one mas) 
would already be very valuable. 
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